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Novel universal correlations in invariant random-matrix models 
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We show that eigenvalue correlations in unitary-invariant ensembles of large random matrices 
adhere to novel universal laws that only depend on a multicriticality of the bulk density of states near 
the soft edge of the spectrum. Our consideration is based on the previously unknown observation 
that genuine density of states and n —point correlation function are completely determined by the 
Dyson’s density analytically continued onto the whole real axis. 
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Random matrices have been introduced in a physical 
context since the works by Wigner JT] and Dyson [2j. Ini¬ 
tially proposed as an effective phenomenological model 
for description of the higher excitations in nuclei j^] , the 
invariant ensembles of large random matrices found nu¬ 
merous applications in very diverse fields of physics such 
as two-dimensional quantum gravity Q], quantum chaos 
H, and mesoscopic physics || . Apparently, this ubiquity 
owes its origin to the very idea of the construction of the 
invariant one-matrix model which only reflects the 
fundamental symmetry (orthogonal, unitary or symplec- 
tic) of underlying physical system/phenomenon but dis¬ 
cards its (irrelevant) microscopic details. Since the sym¬ 
metry constraints follow from the first principles, even a 
rather crude matrix model allows identification of univer¬ 
sal features which persist for a variety of systems with the 
same symmetry. This circumstance emphasizes the im¬ 
portance of the study of universality intrinsic in random 
matrices. 

The simplest invariant random-matrix model is defined 
by the probability density 

P [H] = —— exp {— Tr V [H]} (1) 

of the entries Hj j of the N x N random matrix H, where 
the function V [H] referred to as “confinement poten¬ 
tial” must ensure existence of the partition function Zn, 

1 . 

In the following we restrict our consideration to the 
unitary invariant, U (N), matrix model. Nowadays it 
is widely believed that U (N) invariant ensembles of 
large random matrices with rather strong level confine¬ 
ment may exhibit three different types of locally universal 
eigenlevel correlations which are characterized by the ap¬ 
propriately scaled two-point kernels. 

• Bulk scaling limit is associated with a spectrum range 
where the confinement potential is well-behaved, and 
density of levels can approximately be taken as a con¬ 
stant. It has been proven in Refs. |^| |n| that for rather 
strong confinement potentials |pT|| the two-point kernel 
follows the universal sine law 


A bulk (A S ) 


sin [-7T (s — s')] 
7r (s — s') 


( 2 ) 


Here scaling variable s is measured in the units of the 
mean level spacing: s = e/Ajv- 

• Origin scaling limit deals with that part of the spec¬ 
trum where confinement potential displays logarithmic 
singularity: V (e) —> V (e) — alog|e|. In the vicinity of 
the singularity £ = 0 the two-point kernel takes the uni¬ 
versal Bessel law 1121: 


A-orig (A $ ) 


7T 

~ 2 

( 3 ) 

Here s is scaled by the level spacing near the origin, 
s = e/A N (0). 

• Soft-edge scaling limit, relevant to the tail of eigen¬ 
value support where crossover occurs from a non-zero 
density of states to a vanishing one Jh|, has been only 
investigated for Gaussian unitary ensemble (GUE) [ p~4[ , 
and quite recently for U ( N ) invariant ensembles of large 
random matrices associated with quartic and sextic con¬ 
finement potentials Jn|. It has been found that in the 
soft-edge scaling limit in all these ensembles the two- 
point kernels follow Airy law 

, Ai (s) Ai' (s') — Ai (s') Ai' (s) 

A soft (S, s ) =- - - . (4) 

s — s 

Here s a N 2 / 3 • (e/Djg — 1) with Dn being the endpoint 
of the spectrum. 

Whereas universality in the spectrum bulk and near 
its origin has rigorously been proven for a wide class of 
strong symmetric confinement potentials, the supposed 
universality of the Airy kernel has not been proven. 

Our aim here is to demonstrate that the Airy correla¬ 
tions, Eq. ( p|), being universal for a wide class of matrix 
models Eq. (jjj) , are indeed a particular case of more gen¬ 
eral novel universal correlations which are represented by 
the scaled to— th multicritical two-point kernel 


J Q+ i (tts) J a _ i (ns') - J a _ a (ns) J a+ 1 (ns') 
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K^(s,s') = 


G (a\v*) G' (a'K) • at""* - G (*V) G' (a|i/*) • (s') 5 " 


where the function G is expressed through the Bessel functions as 


G(s| "‘ )= 57t f [ sin (3^) + ( “ 1) 


- 1/2 
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s > 0, 
s < 0, 


( 6 ) 


and parameter v* is determined by the even critical index 
m| m = 0, 2,4, etc. of the matrix model: 

v*=m+Z. (7) 

Note that the critical index m is completely deter¬ 
mined by the type of singularity of the density of states 
near the soft edge 0: vn (e) = (Tr5(e —H)) oc 
(1 -e*/Dl) m+1, \ 

Equations ®) - © together with Eqs. ( p5|) and (j26| ) 
below are the main results of the paper. Although we 
concentrate our attention on the problem of eigenvalue 
correlations near the soft edge, the treatment we present 
here is quite general being relevant to arbitrary spectrum 
range. 

A. Within the orthogonal polynomial technique the 
two-point kernel A/v (e, s') determining the n—point cor¬ 
relation function R n for eigenvalue spectrum of large ran¬ 
dom matrices, R n (e i,...,e„) = det [AW (e;, £j)] ij=1 
can be written through the fictitious “wavefunctions” 
ipn (e) as 


Kn (e, s') = cn 


4>n (s') ipN-1 (e) - i’N (s) tpN-1 {s') 


£’ — £ 


( 8 ) 


Here cjv is the recurrence coefficient entering three-term 
recurrence equation 


sPn — 1 — C-nPn T Sn—lPn— 2 , Pq (e) — 1; Pi (e) — S, 


(9) 


for polynomials P n orthogonal on the whole real axis R: 


I dd (s') Pn (^) Pm (^) — ^nra? 


( 10 ) 


and the “wavefunction” ip n {s) = P n {s) exp {—V (e)}. 
The measure da (e) = exp {— 2V (e)} de is completely de¬ 
termined by symmetric confinement potential 

v ^ = Z& k in) 

fc= i 


To study the eigenvalue correlations in the random ma¬ 
trix ensemble with confinement potential Eq. 0 ) we 
note that a three-term recurrence equation for orthog¬ 
onal polynomials P n (e) can be mapped onto a second 
order differential equation for these orthogonal polyno¬ 
mials and/or corresponding wavefunctions ipn (e)- This 
was already observed for the first time by J. Shohat in 
1930 |T(|. Considerably later Shohat’s idea was devel¬ 
oped by Bonan and Clark ]l7[. The simple and elegant 
method proposed in Refs. @0 turns out to be a very 
general and powerful one for analysis of spectral proper¬ 
ties possessed by large random matrices. 

To map Eq. (|j) onto a second order differential equa¬ 
tion for ip n , we note that the following identity takes 
place: 

dP 

^ =A n (£)P n - 1 -B n {£)P n , (12) 

d£ 

where the functions A n (e) and B n (e) can be found 
from consideration below. Since dP n /d£ is a polyno¬ 
mial of the degree n — 1 it can be represented Jl^] 
through the Fourier expansion in terms of the kernel 
K n (t, e) = J2kZ o p k ( t ) P/t (e) as follows: 

= I da ^ ^r Kn ^ ■ ( 13 ) 

Integrating by part we obtain that 

EST = 2 / da (t) K n (i, s) (A _ Pn ( t ). (14) 

Now, making use of the Christoffel-Darboux theorem |lij| ] 
we conclude that unknown functions A n and B n in Eq. 
( |l2| ) are 

A n (e) = 2c n J da ( t ) —1— —Pn ( t ), (15) 

B n (£) = 2c n J da (t) [> t _l [ ’ Pn (i) Pn-1 (t) . 

( 16 ) 


with d p > 0. The signs of the rest dfc’s can be arbitrary 
but they should lead to an eigenvalue density supported 
on a single connected interval (—Dn, +Dn). 
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At first glance representations Eqs. 0 ) and 0 ) are 
rather useless as far as they involve the same orthogonal 














polynomials which enter Eq. ©• Nevertheless these 
expressions do allow us to get the functions A n and B n 
in closed forms directly related to confinement potential 
and to endpoint of eigenvalue spectrum. Restricting our 
following consideration to large indices n = N 3> 1 we 
reduce Eq. (|J) to the asymptotic form 

ePn = cn (Pn+i + Pn-i) , (17) 

where 

e x P N = Y, ( * ) c nPn+ 2 j-x, A > 0. (18) 

3=0 X J ' 

Substituting V (e) given by Eq. ( llT| ) into Eq. (^b)) yields 

p 2k—1 . 

A n (e) = 2 c N ^2 X! dk£X ~ 1 / da W p n (*) t 2k ~ x ~ l . 

k—1 A=1 ** 

(19) 

Then, taking into account Eq. (|l^) as well as the or¬ 
thogonality of P n we arrive at the expression defined for 
arbitrary e: 


A n (e) = 2c n ^2 dk ^2 ( \ - A ^ 

k =1 A=1 ' 


2k—X2X—2 
C N b 


• ( 20 ) 


It is easy to verify that An (e) can be represented 
through the Dyson’s density 


vd (e) = ~^P 

7 



tdt dV 1 1 — e 2 fD 2 N 
t 2 — e 2 dt y 1 — t 2 jD 2 N 


( 21 ) 


defined on the whole real axis R: 


An (e) 


td'd (e) 
x/l-e 2 /D 2 N 


( 22 ) 


Here e takes arbitrary value that can lie both inside and 
outside of an eigenvalue support. The spectrum endpoint 
Dn = 2cn is the positive root of the integral equation 


M _ 2 [ Dn dV tdt 

TT Jo dt a/ D 2 n - t 2 


(23) 


following from normalization of Dyson’s density. 

Combining Eqs. ©, and using asymptotic 

identity 


Bn = - —, (24) 

Dn de 

which is a consequence of Eqs. (^5[) and (^), it is a 
straightforward step to reach the following remarkable 
asymptotic differential equation: 


<PN- 



nv D (e) 

* 17 rara 


i/j'n + -k 2 v 2 d (e) IpN = 0, 

(25) 


which together with relationship 

* = JZvD-f, 2- ~ ifc**) (26) 

provide a general basis for the study of eigenvalue corre¬ 
lations in arbitrary spectral range [ p~9| . 

An interesting property of these equations is that they 
do not contain confinement potential explicitly, but only 
involve the Dyson’s density vp and spectrum endpoint 
Dn- Moreover, it turns out that the knowledge of 
Dyson’s density (that coincides with real density of states 
only in the spectrum bulk) is sufficient to determine the 
genuine density of states, as well as the n— point corre¬ 
lation function, everywhere. We also note that Eq. (|25| ) 
can be derived in a different way for monotonous con¬ 
finement potentials increasing at least as fast as |e| at 
infinity. This suggests that differential equation Eq. (|25| ) 
should hold generally and not only for confinement po¬ 
tentials having polynomial form Eq. ©■ 

B. Up to this point our derivation was quite general 
without any respect to the soft edge of eigenvalue sup¬ 
port. We now focus our attention on the eigenvalue cor¬ 
relations near the soft edge e = Dn- It is known 0 
that tuning coefficients dk which enter V one can reach a 
situation when the bulk (Dyson’s) density of states will 
possess a singularity of the type: 


vd ( e ) 


1 - 


D 2 
u N J 


1 + 1/2 


TZn 



(27) 


with m = 0, 2,4, etc., and TZn being a well-behaved func¬ 
tion with TZn (1) 7^ 0. [Odd indices to are inconsistent 
with our choice for leading coefficient d p , entering con¬ 
finement potential V (e), be positive in order to keep a 
convergence of integral for partition function Zn in Eq. 
©]• We intend to demonstrate that as long as multicriti¬ 
cality of order to is reached, the eigenvalue correlations in 
the vicinity of the soft edge become universal, and are in¬ 
dependent of the particular potential chosen. The order 
to of the multicriticality is the only parameter which gov¬ 
erns spectral correlations in the soft-edge scaling limit. 

Let us move the spectrum origin to its endpoint Dn, 
making replacement 


e s = Dn 


1 If 2 
+ S ' 2 [ttDnTZn (1) 


i/"*' 


(28) 


that defines the to— th soft-edge scaling limit provided 


oc TV 1 /"*, ft is straightforward 


s < (DnTZn (I)) 171- 
to show from Eqs. (|25|) and (26) that the function 
ipN (s) = ipN (£s — Dn) obeys differential equation 
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$N ( s ) - —- —TpN (s) - s 2( "‘ 1} ^JV (s) = o, (29) 


and that the following relation takes place: 
'IpN—l (s) = VW (s) 


+(-iy 


nD N lZ N (1) 


(30) 

re¬ 


solution to Eq. (|29|) which decreases at s —> +oo (that 
is at far tails of the density of states) is given (up to an 
arbitrary factor X N ) by the function Gfs|i/*), Eq. (!(]). 
The factor Ajv can be found by fitting 0 the density of 
states Km (£ s ,£ s ), Eq. (]|), to the bulk density of states, 
Eq. (|2?j), near the soft edge provided 1 <C s <C -ZV 1 11 /" . 
Then, making use of Eqs. (!), i) and 0 we easily ob¬ 
tain that in the to— th soft-edge scaling limit, Eq. (28), 
the two-point kernel 


A 1rft ( s > s ') = J im K ( £ s< £ s’) -r 1 (31) 

N —»oo as 

is determined by Eq. (|E|). In particular case of to = 
0, that is inherent in random-matrix ensembles with 
monotonous confinement potential, the function G coin¬ 
cides with Airy function, G (s||) = Ai (s), and the Airy 
correlations, Eq. (Q), are recovered. 

It follows from Eqs. (||) and (^) that density of states 
in the same scaling limit 


soft 


(s) = 


d , . t 
Ts G ^ 


-[G(*| 0 ] 


V"* 


(32) 


is also universal. 

The large-|s| behavior of can be deduced from the 
known asymptotic expansions of the Bessel functions: 


(m) 

Aoft 


(s) 


“-(-"A) -’(Jr) 

S in( 3 ^) + (-l)*'*-^ 


5 —> —OO, 

5 —> +00. 

(33) 


Note that the leading order behavior as s —> — oo is con¬ 
sistent with the H" _1 singularity of the bulk density of 
states, Eq. (|27|). 

To conclude, in this paper we presented a general for¬ 
malism for a treatment of eigenlevel correlations in spec¬ 
tra of U ( N ) invariant ensembles of large random ma¬ 
trices with strong level confinement. An important in¬ 
gredient of our analysis is the second order differential 
equation which connects the Dyson’s density with a fic¬ 
titious “wavefunction” ipM which is needed for calcula¬ 
tions of eigenvalue correlations within the framework of 
orthogonal polynomial technique. This consideration is 


relevant to arbitrary energy range. We have applied this 
formalism to examine the eigenlevel correlations near the 
endpoint of single spectrum support. It has been shown 
that in the soft-edge scaling limit there are novel univer¬ 
sal eigenlevel correlations which only depend on the even 
multicritical index of a matrix model. In a particular 
case to = 0, corresponding to monotonous confinement 
potentials, universality of the Airy correlations is recov¬ 
ered. 
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leads to the universal sine kernel in the bulk scaling limit. 
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